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Let  A  be  a  non-negative  n-n  matrix.  In  this  paper  we  the 

growth  of  the  powers  A‘  ,  m  ■  1,2,3,....  .  These  powers  occur  naturally 


in  the  iteration  process 


"^which  is  imjH>rtant  in  applications  and  numerical  techniques.  Roughly 


speaking,  awalyfo  the  asymptotic  behavior  of  each  entry  of  A  ^  •> 

35£Xy'“our>main  result  to  determine  necessary  and  sufficient  conditions 
for  the  convergence  to  the  spectral  radius  of  A  of  certain  ratios 
naturally  associated  with  the  iteration  above. 
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-Significance  and  Explanation 


The  matrix  problem  considered  in  this  paper  arises  when  studying  the 
iteration  process 


where  A  and 
mathemat ical 
other  fields. 


(m*l)  (m) 

x  »  A  x 


(*) 


x*  '  are  non-negative,  '’'his  situation  occurs  often  in 
economics,  population  genetics,  numerical  techniques  and  in 


The  convergence  of  the  above  iterative  procedure  depends  on  the  size  of 
the  largest  eigenvalue  of  A  (the  spectral  radius,  f>(A)).  In  practice  it  is 
laborious  t<>  compute  the  spectral  radius  directly,  and  approximate  methods  for 
estimating  the  spectral  radius  arc  important, 
t'sing  the  notation  of  (•),  let 


.  <»« 1) . 
(x  ) 

,  (x)  «  min  - - - 

i  (*(B,>, 


H  <x)  >»  max 
m  j 


,  (ms  1) . 
(x  ) 

( m ) 

(x  ) 


This  taper  letermines  necessary  and  sufficient  conditions  for  r  (x)  and 
v  (xl  t  tend  t  o  the  spectral  radius  .  (A)  of  A  .  It  is  then  easy  to  compute 

Til 

. (A)  in  practice.  The  above  problem  is  completely  analyzed  by  determining 
the  a.symptot  i c  behaviour  of  each  entry  Of  A™  . 


For  any  non-negative  x  denote  r(x)  «*  min(Ax)./x  and  p(x)  *  max(Ax)  /x 

l  1  i  i  1 


We  determine  necessary  and  sufficient  conditions  in  terms  of  the  reduced  graph 

of  A  such  that 

lim  rtA^x)  ■  *(aWx>  *  o(A)  . 
m  •**  m 

This  Is  important  for  numerical  procedure  for  calculating  o(A)  -  the  spectral 
radius  of  A  . 


T  re  s;  'ir  cib  i  1  i  ty  for  the  wording  and  views  expressed  in  this  descriptive 
ummar  ’  lies  with  MFC,  and  not  with  the  authors  of  this  rej>ort. 
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1 .  introduct  ton. 

Let  A  be  a  non -no)  at  ive  n*n  matrix.  In  the  iteration  proce** 

U.l>  A,  .‘0)>0. 

which  1*  important  in  app  1  icat  ion*  and  numerical  techniques,  the  |ow»r» 
occur  naturally.  In  thia  paper  we  »tudy  the  growth  of  these  power*,  when  A  i*  irreduc¬ 
ible  or  *tocha»tlc,  the  behavior  of  A*  i*  well  studied,  e.g.  tiantmacher  [7 ,  Ch .  IJ,  iS-7)  , 
Varga  IP,  pp.  12-14).  In  these  case*,  the  eleeentary  divisor*  belonging  to  the  spectral 
radius  p(A)  of  A  are  linear.  We  deal  here  with  the  general  non-negative  caae,  when 
the  elementary  divisors  belonging  to  f (A)  may  have  degrees  greater  than  1  .  At  the 
ost  of  Ignoring  nllpotent  A  ,  where  the  problnm  is  trivial,  we  auvr  that  p(A)  -  1. 

For  a  complex  n*n  matrix  A  ,  with  p(A)  •  1,  there  is  a  least  integer  k  for 
-k  m 


However,  even  in  the  simple  caae  of  an 
-k  m.  *  ....  . .  -k  _m 


which  a  ~  A~  is  bounded,  m  -  1,2,1,... 

imprlmitlve,  irreducible  non-negative  A  ,  1  im  Im  *  A*;  and,  a  fortiori,  11m  »  A*  , 

m  —  m  — 

So  not  in  general  exist.  To  obtain  precise  results  fat  general  non-negat  ive  k  with 
r(A)  •  l  It  la  thus  necessary  to  introduce  sosse  smoothing.  For  example,  in  III)  Pothblun 
considered  Cesaro  means  of  powers  of  A  .  In  this  f>ap>er  we  study  the  growth  of 
a. 2)  It1*'  -  A*  <1  ♦  ...  .  A<,-*>.  m  .  1.2,... 

where  q  is  a  certain  positive  integer. 

After  s.**e  preliminaries  in  <2,  we  use  elementary  analytic  methods  in  I)  to  prove 
a  theorem  on  the  grtmeth  of  H  .  As  corollary,  we  obtain  a  known  theorem  on  the  index 
of  the  eigenvalue  1  of  A  ,  cf .  Schae'er  Ilk,  Ch .  1,  7S».  2.7|.  We  also  give  a  local 
fora  of  the  theorem,  that  is,  we  show  that  for  1  •  i,  j  •  n  there  exist  integers 


k  -  k(i,J)  and  q  •  qd.jj  0  such  that  the  el  resent  b 


(m) 

i) 


of  the  matrix  given  by  (1.2) 


sat isf ies 

(1.)) 


lim 


-A  _<■> 

*  bil  >0- 


The  analytic  results  of  *1  motivate  the  investigations  in  the  rest  of  the  paper. 


on  sored  by  the  ttnited  States  Army  under  Contract  Mo.  r  >AA^29-7$-C-0024 .  This  material 
Is  based  upon  work  supported  by  the  national  Science  Foundation  under  Crant  Mo.  MCS7B-r  1  oft? . 


The  n.*in  thrust  of  the  paper  is  the  use  o!  the  qraph  structure  of  tlie  Matrix  A  to 
decrease  the  intoqet  q(l,J)  and  to  determine  the  mteqer  k(i,J)  In  (1.1).  The  req¬ 
uisite  jiaph  theoretic  .01  cjt#  are  developed  in  14  and  in  i‘>  we  state  our  Main  result. 
Theorem  (S.l  ).  As  a  corollary,  we  ottain  a  striking  theorem  on  the  index  of  1  due  to 
Rothfclua  113). 

tn  Sti  we  apply  our  results  to  the  Iteration  process  (1.1).  For  x  0  ,  x  x 
denote 

(1.4.1)  r(x)  «  eupfuiux  ^  Ax'  , 

(1.4.11)  X(x)  •  Inf (pipx  ^  Ax ! . 

SI  g# 

In  Theorem  (6.8)  we  find  necessary  and  sufficient  conditions  for  r  (A  x!  and  ft(A  x)  to 
converge  to  the  spectral  radius  of  A  .  In  )7  ,  we  show  that  a  theorem  due  to  D.  It. 
Carlson  (1)  on  the  existence  of  non-neqat l ve  solutions  y  for  (I  -  A)y  *  x  ,  x  •  0  is  a 


consequence  of  our  main  results  and  we  extend  the  theorem. 


be  a  imiimac*  of  non-negat  1  vo  and  k  _  0  be  an  integer 


i  a  bounded 


will  denote  that  v  (■)/«• 


will  denote  that 


will  denote  that  ll»  exists  and  la  jmsltive 


he  stove  notation*  will  alto  be  ueed  for  k 


will  each  Indicate  that  there  exist*  r 


the  above  notation*  will  Bean 


P  <  1  ,  *uch  that  <H*Oc> 


( 1) .  In  case  that  k 


that  *(0)  *  1  for  all  sufficiently  large  »  .  (Thu*  «t(») 


(v)  The  notation  Mb)  «•  »  will  be  u*ed  'or  a  sequence  of  non-negative  matrices 


tf  mil  ate  the  relation  hold*  for  each  element 


CanMnator  ial  knuit 


Let  t  2.  0  and  t  >  0  be  integer*.  Then 


<2.21  t  *  L  1  I  ...  x 

V’  •  *pt“r 

where  'he  vaMtuw  is  taken  over  all  non-negative  integer*  p 


It  1*  well  known  that 


The  six:  leu*  way  to  prove  thia  equality  is  by  considering  the  coefficient  of  x  of  both 


•ides  of  the  identity 


which  i*  derived  from  (1-x) 


for  a  purely  combinatorial  proof 


see  for  ex**|  1  *  Brualdi  12.  p.  J7)  .  for  t  -  0  the  above  formula  implies 


We  Shall  also  need  »rw>e  result*  on  the  convergence  of  eerie* 


t 


(2.4)  I,*a»a ;  Let  k  ^  1  and  let  b^  0,  p  -  0,1,2,...  be  a  sequenco  such  that 


(2.5) 

where  q  »  0  .  Then 

(2.6) 


-<k-i)  , ,  , 

p  lia  (b  ♦  . . .  ♦  b  ,)”V, 

p  -  -  P  P**-1 


a  k  lia  b  •  . 

b--  p-1  P  *» 


•k  k-1 

Proof:  Eleeentary.  klternat i vely ,  check  that  c  •  »  k  p  satisfies  the 

a,p 

assumptions  of  Marty  (0,  Theorem  2,  p.  43]. 

(2.7)  :<rmA  i  Suppose  (2.5)  holds.  If  lia  a  -  u  then 


(2.8) 


-“it*  [  -  _ 

■  **  p-1  P  "'P  *“» 


Proof  :  According  to  Hardy  (8,  Theorem  16,  p.  641 


1  *b 

P  "*P 


(2.9) 


since 


t  IB 


J=1 


l  > 

*•  _  n 


i-i 


o  < 


I  b 
p^i  p 


b  ♦  ...  «  b  .  ,  (k-1) 

a _ ’ _ a»g-l  2  v  a 


i  b 

1-1  P 


v(2kq)*1  a* 


ar>d  the  last  expression  tend*  to  0  .  If  w*  apply  (2.6)  to  (2.9)  we  obtain  (2.8). 


By  3-  ,  reap,  ff  ,  w  denote  the  real,  reap.  complex  field,  and  by  D-^  the  non- 

negat  lve  numlfts.  The  aet  uf  real.  reap.  complex,  non-neqative  r*n  matrices  will  be 

denoted  by  D*'  ,  reap.  Cf'  ,  1F^  '  .  Me  alao  write  A  ^  0  for  A  c  B*'  (A  i»  non- 

neqative)  and  A  ■»  0  when  A  la  poaitive  (a^  »  0  ,  i  •  l,...,r,  )  •  l,...,n). 

Let  A  •  rnT!.  By  apec  A  we  denote  the  aet  of  ciqenvaluea  of  A  .  Suppoae  that 

a  pec  A  •  <1  . I  >  ,  where  the  1  are  pairwiae  distinct.  It  ia  known,  cf .  (Jantmscher 

(7,  *h.  S,  1 3 )  ,  that  there  ealat  non-neqative  inteqera  Pj....,pjr  and  unique  matrices 

,(  i.)  (  pnn  ^  g  .  . , . ,p  ,  i  *  l,...,r  which  are  linearly  independent  auch  that  for 

9 

each  polynomial  fit) 


The  Z*  ''  are  po  1  ynern lain  in  A  ,  p^  ♦  1  is  the  size  of  a  larqest  Jordan-Mock 

tap) 

helonqlnq  to  i  .  The  coliaana  of  Z  1  are  eigenvectors  of  A  corresponding  to  the 
3 

«•> 

eigenvalue  V  ,  the  rank  of  7.  la  equal  to  the  nisaber  Jordan  blocks  of  sice 

3 

p  ♦  l  corresponding  to  i  (The  aim;  lest  way  to  obtain  (3.1)  ia  by  assuming  that  A 

3  3 

is  in  Jordan  form).  Aa  usual  we  define 

indeed  )  •  p  *  1. 
i  a 

That  la  p  ♦  l  la  the  multiplicity  of  i  in  the  minimal  polynomial  of  A  .  We  shall 
3  3 

also  use  a  local  iced  index.  Tor  1  ^  1,J  ^  n  we  put 

index j  ^ V^)  *  1  *  maxleiaj’  *  0  ,  f  •  0,...,p)  , 

where  index  O  )  •  0  if  t."*'*  ■  0  ,  t  *  0 . .  .  If  At  *nn  and  m  is  any 

» J  3  1 )  3 

integer  we  shall  denote  the  elements  of  a"  by  a|**,  1  i,)  m  . 

Let  A  t  ff^'r  Me  assume  throughout  the  normalisation  c(A)  •  1  .  It  ia  well- 

known  (Frobenius  |*>1  .  "-antmarher  |7,  ch .  13),  Berman-Fl emmons  |1,  Ch .  I) )  that  if  >  is 

an  eigenvalue  of  A  and  l|  *  1  ,  then  A  la  a  root  of  1  .  Hence  there  ia  a  positive 

integer  q  such  that  >**  •  1  ,  for  all  )  t  spec  A  ,  !»'  »  1  .  The  smallest  auch  integer 

q  will  be  called  the  r<r  iod  of  A.  If  q  »  1  ,  A  will  be  called  aperiodic.  For  ah 

Irreducible  and  aperiodic  naT rix  A  •»  0  the  Frobenius  theorem  and  the  formula  13.1)  imply 


ll*  A 


,(10) 


>  0 


where  .  j  »  1,  set  for  example  9cm*n-P!ra»w«  11,  Tt*».  4.11.  Theorem  (1.4)  extends  the 
above  equality  In  a  local  way.  hart  (l)  of  the  theorem  is  an  extension  of  the  known  in¬ 
equality  apparently  due  to  Schaefer  114,  p.  264,  Th* .  2.4) 

0.2)  index  (1)  <_  index  (l)  ,  If  |l|  •  1  , 

for  non-noqat ive  eatrices.  see  also  Schaefer  [IS,  Ch.  I,  Ttm.  2.71,  Beraan-Plramnns 

(1,  Thm.  1.21.  This  result  and  Part  (i)  of  Theorem  (1.4)  could  easily  be  deduced  from 

the  classical  pringsheia  theorem  functions,  e.q.  Tltc'marsh  [16,  p.  2141.  The  use  of  the 

Pringshctm  theoree  in  analyzing  the  s|>ectral  properties  of  non-ne<jat  ive  matrices  can  )>e 

traced  back  to  atrowakl  [10),  see  also  Farlin  (9)  and  Schaefer  [14.  App . 1  for  the 

infinite  dimensional  case.  See  rnedland  (S)  for  a  detailed  analysis  of  the  Pringsheim 

theoree  for  rational  functions  which  has  certain  analogs  of  the  Frobenius  theoree.  For 

sake  of  .xaipelteness  we  brinq  a  short  and  elementary  independent  proof  of  Theorem  (1.4). 

To  do  so  we  need  an  easy  leessa  which  probably  is  known. 

0.1)  heama.  bat  ,t  .  >  -  1 . r  be  cosip  lex  numbers,  where 

-  s  a  r 

the  >  are  pairwise  distinct.  1*  1  in  (  /  >  z  )  exists,  then  z  -  0  if 

a  *•.  a  o  <1 

»  •  •  '*•1 

ii  >  i  ,  i  *  i  . 

3  —  <1 


Proo 


fi  Since  lisi  »*  exists  for  !i  1  or  1 


we  sia\  assume  that  i  >  1 ,  d  1,  a  •  Put 

r  —  o 

u*  —  (u  , . . . ,  u  ,).  where  u  *  z  ,  let  1  * 

•  «.r-l  »  aii  « 

V  •  (v  )r  -  C,r  he  the  Vanderaond  matrix  given  by  w  * 
aB  1  of 

U  ■  W(  t  , 


The  assissption  of  the  !iew  implies  that  1  im  u 

■  •  • 

Its  1  i  •  Us  v  u 


(a) 


exists 


•  1  ,  without  loss  of  generality 
*  •  (Zj . zf)  <  f  ,  and 

dlagllj,...,!^  .  Jrr  and  let 
1*1 

a  f  7h«»n 


Since  V  is  non-singular. 


and  so  z  •  0  . 


0.4)  Thcort 


Let  A  •  D-^  ’  where  c  (A)  »  1  .  Let  1  ‘  i,  )  ^  n. 


(I)  If  i  .  spec  A,  1  V  |  •  1  ,  then  indexed)  ^  Indexed)  . 

(II)  Let  q  be  a  positive  integer  such  that  i'*  •  1  if  X  ,  spec  A,  |  A  |  •  1  and  let 


k  ♦  1  •  indexed).  Let 


.<■) 


A*(I  ♦  ...  ♦  A«'1) 


Then  b 


(•) 

1} 


In  ■'articular,  a^**  4  o(aN  ,  if  k  0 


Proof i  (i)  Let  }  be  the  eigenvalue*  with  |i  ■  •  I  ,  n  • 

—  ■-  I  r  a 

where  the  >  are  pairwixe  dixtlnct.  Let 


d  ♦  1  •  aax( index  .  (1  ):  a  •  l,...,r  ). 
ij  a 

If  d  •  -1  then  there  ia  nothinq  to  prove.  So  asstsse  that  d  _>  0 

Suppose  that  r^  r*  '‘'eO  for  a  ■  where  1  <_  s  <_  r  .  It  follows  isssediately 

f rcei  (J.l)  that 


Hence,  by  L 


.  „«<  y  i*-d  ,  ,  ♦  0(.d)  . 

l)  a-l  °  a 

O.  J)  ,  aj*1  4  o(od)  . 

Let  q  be  a  positive  integer  such  that  i*!  *  1 

9 

4  (T)  -  t"(l  ♦  t  ♦  ...  ♦  tq'l>. 


1 , . . . ,  s 


Define 


If  w*  take  the  d-th  derivative  o'  4  (t),  we  obtain 


<i  , ,  ,  <* 

*  4  .<  t)  ♦  o(a  )  , 

■  -d 


for  any  fixed  t.  -  ’  1  ,  and  also  -  0  for  |>J«1  i  fl,  1  <  a  <  s.  put 


.  (■) 


B  -  By  0.0  and  the  equality  above  we  have 


0.5) 


U 


(  f  V  A(*  )t  )  ♦  o(p‘)  . 
*•,  ■-<!  i  l 

9*1 


Vow  suppose  that  indexed)  <  d  ♦  1.  Then  O.S)  Implies  that  bj"*  -  o<»'1)  .  But 


(>» 


(■) 


,  “  •.  .  —  0  Bnd  this  is  a  contradict  Ion.  Thus  d«k  and  this  proves  (i) 

(1‘)  Suppose  that  i  •  1  .  By  0.5)  and  the  preceedlnq  arguaent  we  obtain 

.(•)  k  .  k. 

bl)  *  •  **l  4  °<B  ’• 

where  i,  •  ►  0  .  This  prove*  (11).  • 

we  now  state  a  qlobal  version  o'  Theoreo  0.4) .  oart  (11)  which  follows  iamediately 


frcai  Theoresi  0.4). 


0.6)  Theorsa. 


Let  A  .  D- 


where 


p  (A)  -  1  . 


Let  q  be  a  poiltive  lnteqer  such 


that  »*• 


Then 


1  if  V  «  spec  A,  !  X  |  ■  1  and  k  *  1  •  index  (D  -  index  ID.  Let 


„(a)  .q-l. 


ii.  a*k  .'■»  -  r  . 


where  r  ^  0  and  y  is  not  Identically  tero. 


4.  ' »  .-%!  t.  theoretical  ept  s  . 

let  A  '  V  ''  and  let  .(A)  *  1  .  We  may  assume,  without  loss  of  generality,  that 
after  simultaneous  permutation*  of  rows  and  columns,  A  is  in  the  rrobenius  16)  normal 
f  ■  >  rts  which  can  1*  fimnd  in  many  references  e.g.  Gantmachn  |7,  Vol .  II.  p.  7  5)  .  Thus 


(4.1) 


11  "12 


22 


'lv 


where  the  diagonal  blocks  A  ,  a  •  l,....  v,  are  irreducible  and  all  subdiagonal 

33 

blocks  are  0  .  (The  1  •  1  matrix  0  is  considered  to  be  irreducible). 

bet  A  he  in  Frobenius  normal  form  (4.1).  Then  the  (reduced)  g r ap  h  (1(A)  of  A 

is  a  subset  of  <v>  *  <v>  ,  where  (v>  •  (l,...,v)  and  C(A)  •  <(3,6)  «  <v>»  (v>  :A  4  0). 

as 

(Observe  that  many  authors  w<sild  call  (1(A)  the  arcset  o*  the  graph  «.  >,C(A)),  but  we 
have  no  nee<)  to  mention  the  vertex  set  (v)  explicitly). 


If 

la, 8)  «  C (A)  we  call 

(a.0> 

an  arc  of  C (A) . 

If  (3,6) 

is  an  arc  of  CIA)  , 

then  a  ^ 

6  and  (3,3)  <  CIA)  , 

1  ^  a 

< 

V  8 

unless  A 

33 

is  the  1 

*  1  matrix  0  .  Thus 

we  define 

a  (simple)  iath  f ror. 

a  to 

0 

in 

C  (A)  to  be  a 

segoence 

s  »  (aQ . O  ,  where 

either  s 

^  1  ,  1  3  •  0Q  <  ... 

<  a  • 

• 

0 

L  v 

and  I cst _j , 

ij)  *  C (A) 

• 

O 

» 

• 

A— 

or  s  •  0 

and  a  •  o„  •  8  and 

0 

(o,o) 

r 

C(A) 

The  suj^xj 

it  o£  * 

is  the  set  »upp  *  • 

t  (l,...,v).  We  always  assume  that  the  1^,  i  •  0 . s  liave  been  listed  in 

st  r  1  t  l  v  ascending  order  . 

If  l  <  3  *  v  ,  then  we  call  1  a  singular  vertex  (of  0(A))  If  r  (A  )  «  1  . 
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(This  terminology  is  •-(«,<» i *>n*  with  that  of  ^ lchman-Schnelder  111)),  bet  1  3  ^  r  ^  v. 

Tor  any  path  -  frrss  n  to  8  let  k(*)  ♦  1  be  the  nussber  of  singular  a  in  the  sup¬ 
port  of  s  .  (Thus  note  each  distinct  i  is  counted  only  once  in  k(e)  *  1).  bet 

3  <3  s  ...  %  ,  where  k  »  k(e)  ,  be  sll  singular  vertices  in  supp  »  .  If  there 

*1  ’k 

is  no  path  frrss  ,3  to  ?  in  '.(»)  we  put  k(s)  .  Noxt  we  set 

(4.2)  k(i,f)  *  max!k(»)!  «  is  a  nath  fr<»s  1  to  p  in  GtA)  I. 


-g. 


We  ahall  call  *1  !,.•*>  the  singular  'll  at  «t»r  frusi  i  to  #  .  If  (l,i)  1»  a  position 

in  A  and  (}.))  a  position  in  A  then  we  ahall  also  call  k|».)1  the  singular 
11  rf 

distance  frt«i  1  to  ]  (note  our  u»e  *>f  square  brackets). 

A  lath  •  fj  m  ■  to  :  will  >*•  ailed  a  naaiaal  ;  at  h  If  the  inaler  of  singular 

vertices  In  the  supjwrt  of  *  la  kl  t,R)  •  1.  U>t  I  ^  j,  I  ^  .  bet  P(a,b)  be  the 

aet  of  saaaiBal  jatha  frtwa  a  to  t  .  hor  each  »  «  Pla,4)  let  qlt)  be  the  q.C.d.  of 
t'orlod#  of  A  with  i  <  eui  f'  v  4 ml  singular  (vti.  ,  (A  J  *  II. 

»»  rf 

than  we  define 

(4.J)  q(i, j)  •  l.c.a  q(*)t»  >  Pfa.tl). 

We  shall  all  q  (  i .  r  I  the  1  --  a  1  i  er  ic«l  of  (a,:).  If  k(i.f)  0  t  lien  <|  (  i ,  f  1  •  1 . 

Also  If  (l.l)  la  a  position  in  A  and  (),})  la  a  position  In  A 

At!  P  P 

pot  q(*,f)  •  q ( 1 .  )) .  the  local  period  of 


t  hen  wc  ahall 


5.  The  main  losuito. 


Let  a  l>t  ,  whole  i  (A)  •  1  .  lw  in  rtobtnlui  norwl  for*  <4.1).  It  follows 
ftijm  the  rtrroti'frohtniM  theory  for  non- negative  matrices,  e.g.,  Gantmacher  (7,  ch.  11) 
that  t  here  is  a  dlaq  nal  nat»i*  t  with  positive  diag.wial  clcesents  so  t  )at  ,  upon  re)  lacing 
A  by  XAX 

(5.1)  A  •  p(A  ) A*  , 

ta  3!t  aa 

whore  A '  is  a  stochastic  mat  ns, 

Ml 

15.2)  IA  t  <o,  1  <  a  <  |  ‘  v  ( 

IF  •  ™  —  — 

where  1  »  i  Ml '  (A  )t  t  (A  )  *  1  ,  i  •  l,...,v).  Mere  (I  is  the  i  -operator 

»  li  •  • 

no  ns . 

n 

1**^  •  maxi  «  !r  1  •  1 . r)  for  X  c  * 

l*1  5  ,  , 

The  diagonal  »»tni  X  can  le  constructed  as  follows.  bet  u  lie  a  positive  vector 

#At  laf  y  ing  A  u*  '  •  ?(A  )u!l  .  Denote  by  X  a  diagonal  matrix,  whose  diagonal 
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entries  are  the  elements  of  u*’*.  Then  X  is  of  the  form  dlag  ‘  X  ,  <  X  , . . . .  i  *X  )  for 

12  v 

« rate  small  enough  positive  i  .  In  our  subsequent  proofs  we  may  assume  that  A  has  been 
normalised  as  shove. 


Ut  t 

be  a 

path  in  CIA)  .  Denote  by  a  *  1  the  cardinality  of  supp  *  . 

(S.J.l) 

supp  . e  }  , 

*  i  *0  *  ^3  *  ...  *  1,  i  v. 

We  def ina  the 

.V  ’ 

matrix  A(»t  by 

Aii<*)  "  \$k 

i  •  0, . . . ,s 

(5.1. ID 

‘•1#I  Vlel 

i  •  0, . . . ,s-l 

Aj  ^  (s)  •  0,  otherwise 

1.)  *  o . . 

(5.1.111) 

Al»)  -  (A  (*) ) *  . 

ij  0 

Thus  A(»>  is  in  Trobeniu*  normal  form  and 

has  *•!  irreducible  diagonal  blocks 

“A  ,  i  •  . . .  To  avoid  ambiguity,  we  write  A(*)**'  for  the  (1,J) 

11  ®l*l  *’ 

block  exponent  of  A(*l  .  i,J  •  0 . S. 

we  now  prove  a  aequer.ee  of  '.«*>«<  for  the  rath  matrix  Alt)  of  a  given  path. 


(5.4! 


i.  Let  A  ,  I*  where  p(A)  *  l.  Let  1  a,  6  ^  v  end  »  be  *  path  in 


G(A)  ti<m  4  to  t  .  Put  A  •  A(»),  where  A(»)  •  l  is  the  masher  of  unquiet  vertices 
in  supp  »  .  If  A(*>  is  the  path  aatrii  given  by  (5.1),  then  fA(»)<*,f-  •  0(*k)  . 


(5.5) 


Proo£:  We  note  thet 

«■)  T 


At*) 


Os 


♦P. 


\>,Aol,’,Au(',"A.-l).,',\>  • 


So 


lA(»)_  I  *0 

0»  *  - 


p  •  p 

t  A  (»)l  °  ...  IA  (»)l  * 

oo  •  •»  * 


p0 . *V 


Suppose  firet  thet  »  does  not  contain  singular  vertlcee.  l.e.  A  •  -1.  Then 

e  »•! 


|A(e)‘-’.  .  c* 

Os  •  — 


l 


1P°...1P- 


P  ♦ . . e*P 

0  • 

where  f*  is  given  by  (3.J).  Ae  r^”*  ^  a*  we  laaedistely  deduce 
lla  t"*A(»)*>  -  0  ,  for  any  t,o  <  t  *  1  . 


Suppose  now  that  A  ^  0  .  Than 


IA(«)!"’»  <  0*  l  \  °...l 

0a  •  —  4 


S  Vi 

o  . . .  o 


q  ♦  •  •  •♦<!  *•*• 

0 


**  r  <  i 

ueO  q0» . . .eqk»u 


iq°...l,k) ( 


l 


■a.1 


»-• 


*  rU 

9  I  rkel  i-A 


Hence 


IA1*’  (e)  I  • 


u»0 


OS  -  —  A»  1 


‘  L  r.-A 


V 

0  )  a 


N^O 


The  last  series  converges  by  the  retio  test  and  j  *  This  establishes  the  leap*.  • 

(5.6)  Let  the  assuapt  ions  of  tesasa  5.4  hold.  furthensore  that  A  ^  0  , 

l.e.  the  support  of  e  contains  singular  vertices.  Then,  for  sufficiently  large  a 
2(s*l) (n-l) 


(5.1) 


T 

)-0 


A(*l 


*-»»  > 


oc 


where  is  a  positive  natrtx. 


*12* 


t 


Proof  :  Let 


we  have  b  (a)  >  0 


Since  A  (»)  is  irreducible,  end  its  distension  does  not  exceed  n 


Nlelandt  118),  bernan-Pl 


s-l.s  ••  M  s* 

t»e  the  natrix  all  of  whose  entries  equal  1  and  whose 


linens  ion  is  that  of  A  (»).  Clearly  B  <»)  ^  CAE 


s-1,  and  hence  for  scale  c  »  0 


where  c.  »  0 


In  the  inequality  (5.">  we  nay  restrict  the  s>ai  on  the  rl<jht  hand  side  by  lettlnq  p  *  0 


*  i  be  the  subscript*  of  B  which  are  sinqulsr 


it  follows  that 


vert  ices  and  put  A 


are  tsat  rices  all  of  whose  entries  are  1 


k .  It  follows  that 


for  sufficiently 


lerqe  • 


write  faqain  without  confusion)  index,  ..(11  *  nax  lndex  (11 
^  (a,fl  i.) 


(%.<»)  u 


tot  the  asauaptiona  of  L 


(5.4)  hold.  tot  q<*>  he  the  g.c.d.  of 


per  loda  of  A  for  aingular  i  aupp  ’  .  tet 

Bit)1*'  •  »(*)"(!  •  M»)  «  ...  ♦  Alt)1'  1 )  , 


4  •  q(  *) .  Then 

<11  In  At»),  index  (1)  •  k(»)  ♦  l  > 

O.a 

Cli)  1  is  a)  !*'  >  0  . 

k--  0'* 

Proof  i  let.  •  ♦  1  •  index  ^(1)  in  A(»).  By  Theorao  (1.4)  there  la  an  integer 
q*  ,  l  _  q*  •  n  ,  nmh  that,  for 

C<»)<-’  -  A(a)<B>II  .  A(»)  •  ...  ♦  A(»)q*'1), 
ju*  *  0.  If  k^O. 

But  .  k(t)  -  k  intradicta  :.nma  (V4),  and  •  •  k  contradicts  l.nema  (S.fc).  Hence 
«  •  k.  Thla  prove#  (1)  . 

flow  let  q  q(»).  If  1  ia  in  the  apectru*  of  A(»),  *  •  1  and,  in  A<»), 

index  .{•)  -  •  (»'  •  l,  then  «n»at  l»  an  eigenvalue  of  every  A  ,  for  aingular 
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-t  «upp(»).  Hence  *  -  l .  q  -  ql-).  He  imediately  obtain  (ii)  free#  Theoreaa  (1.4)  and 
l-vma  < J.6) .  • 

Vo  now  atate  cor  nain  remit. 

(k.l  )  The  r«e»  ■  let  A  f*  nor-irro  n-n  eatiix  normal  ited  by  the  condition  o(M  •  1. 

Aaa.jae  1  _  l,  '  _n.  let  k  •  kli,))  be  the  alnqular  diatance  f roe  i  to  )  and 

q  *  qli.jl  be  the  local  period  of  (I,)).  Put  B1"'  -  A*(l  •  A  »  ...***  *),  Then 


Proof  Ax  uaual.  w»  aaavjee  that  A  ia  In  the  Probeniua  fore  (4.1). 
(i.i)  i*  a  poaition  in  A  and  (),})  a  poaition  in  A..  .  Denote  by 
*»«»?  of  «ll  paths  oonn* -•  imj  ■»  to  ?  .  Th*n  v«*  obviously  have 

.<■>  r  ...(■) 


s'v  -  i 

'*  a.flia.t)  C»‘* 


Suppoae  that 
It  (a.  8)  the 


So 


l 

*'■(•>, S> 


B(t> 


(■) 

0a(») 


Aaaveie  firat  that  k  *  k(»)  •  ••  then  clearly 

lama  'S.4)  iapliaa  that  each  A(*)**'  .  *a  *. 

0*(*) 


So 


(■) 

•  k  x 
at 

- 

at 


«  0. 

-1 

■ 


If  k  •  -1  ^  k(»)  then 

and  again  A1*'  •  B,b!  . 

at  at 


-14- 


However,  If  k  »  kl»),  then  according  to 


a*  i}(»)  divides  q(.i,|5)  *  :j  I  i , )  I  .  By  the  definition  of  kin, 6)  there  exista 

.  Ml.  Ml  _l  II.I  blal  _  L  I  -  .  I  Ml  11.  ”  ( X  )  Ml 


■nder  the  condition*  of  Theorem  (A. 10) 


If  k  •  -l  ,  then 


I’roof  i  for  k  •  0 


by  Theorem  < 5 .  10)  the  non*nogetlve  aeries  above  converge*.  The  aasuaptlon  k  •  -1 


labile*  that  at  leaat  one  tern  la  positive,  finally  if  k 


and  the  result  follows 


.'tapering  Theorems  J.4  and  S.la  ve  first  deduce  a  local  version  of  Rothblun'a 


equality  and  then  the  equality  itself 


then 


where 


Index  (1) 


(i.lll  Corollary.  (hothblum  till)  let  A  P 


Then 


index  (1) 


index (l) 


Evidently  our  result*  <5.10)-<S.D)  are  easily  modified  to  apply  to  all  non-negative 


6.  Convergent  itcretive  methods  for  the  spectral  radius  of  a  non-negat  1  uc  muu, 

A  «  *»""  and  »«•».  that  0 <X)  >  0  .  Le*  r(x)  and  Ml)  t-e  defined  ai 
(X.4>,  Clearly  0  <_  r(x)  <  B(x)  <,  •  •.  It  is  obvious  that 

r (x)  <  r (Ax)  <  r (Ax)  «  R(x)  . 


So  tha  sequence  r (A  x) ,  »  «  0,1, 


It  an  increasing  sequence  bounded  above  by  Six) 


and  the  sequence  R(A  x) ,  a  •  0,1 .  it  a  decreasing  sequence  bounded  below  by  r(x). 

In  [41  Col  1st*  observed  that  for  A  >  Hnn  and  x  »  0 

♦ 

(6.1)  r(x)  <_  »- C A)  *_  Mi)  , 

and  when  A  is  irreducible,  this  inequality  is  valid  for  all  x  ^  0,  x  ^  0,  see 

Mielandt  (1B|,  Varga  (17,  p.  121.  Thus  tha  question  arises  when  for  A  ^  0  and  x  0, 


lie  r(A*x)  e  c (A)  e  1  IB  R(A*X)  . 


Wlelamlt'a  (IB)  characterisation  of  o<A)  for  irreducible  A  easily  implies  that 
(6.2)  holds  for  primitive  A  and  all  x  t  x  ^  0,  x >  0  cf.  Varga  (17,  p.  34).  This 
result  follows  Ire**  the  fact  that 

1  la  c(A)**A*  -  t  »  0  , 

when  A  is  primitive,  where  Z  •  uv* ,  v  »  0  Au  «  c(A)u,  v  >  0  ,  v*A  •  pfAlv*. 
v*u  ”  1 .  If  A  is  irreducible  but  imprlmitive  then  (6.2)  does  not  hold  unless  x  is 
orthogonal  on  all  eigenvectors  of  A*  corresponding  to  >  such  that  jij  •  o(A)  and 
>  4  r(A>.  v*  shall  show  that  this  condition  can  be  cut  in  equivalent  forms.  If  A  lr 
Irreducible  and  of  period  q  ,  then  by  eloultaneous  oermutat ions  of  rows  and  coliasns  v* 
now  put  A  into  the  form 


where  the  diagonal  blocks  are  square ,  rrobemus  (6),  Gantmacher  (7,  Vol  II,  p.  621, 

Herman  -Fl<«>on>  (1,  Ch.  2,  Tho .  2.20), 


(6.4) 


Let  A  be  an  Irreducible  non-negative  matrix  of  period  q  In  font  (6.1), 


and  suppose  that  p(A)  *  l.  Let  v*A  •  v' ,  Au  ■  u,  where  v  >  0  ,  u  •  0  ,  v*u  *  1, 


A*y5  -  J,.',  )•  l,...,q-i,  w  •  e*  l/<*  .  Let  0  *  s  <  be  partitioned  conformally 


with  A  ,  s'  •  fa* jj  )  .  Then  the  following  are  equivalent 


(i) 

(11) 

(ill) 

(iw) 

(w) 


lim  A  a  •  (v  a)u  , 

■  •  « 

lim  As  eaiata  , 

m  ♦  • 


»V-o.  , 


t 

V  s 
(1)  (1) 


1 . q-l 

t 


"  (q)  (q)  * 

Urn  MA%)  "lim  r(ABs)  »  1, 


where  v 


(v 


(l)' 


,v(^j)  has  been  partitioned  conformally  with  A  . 


■roof :  We  first  derive  a  formula  for  As,  m  *  1,2,...  .  Let  w  be  a  primitive 


q-th  root  of  unity.  It  is  well  known  that  the  eigenvalues  of  A  on  the  unit  circle  are 


a-1 


i  *  »'  ,i*l . q  and  that  each  *  is  a  simple  rero  of  the  characteristic  poly- 

a  a 


nrssial.  It  follcars,  in  the  notation  of  11,  that  p  *0,  a  •  1 , . . .  ,q  and  that 

a 


,  (aOl  a-1  tl-a 

7.  •  0  uv  a  ,  a  *  1 ....  ,q 


y  »  0 


(1-a) 


a  »  1 1 


.q-1  . 


where 


It 


wl 


22 


and  I  is  an  identity  matrix  of  the  same  order  of  A 

M  aa 


•  1 , . . . ,g  . 


Hence  by  (1.1) 


Y 

a«0 


11  t-  a  ... 

n  uv  o  *od) 


and  so 


-17- 


J 


Ax*  )  •  4  ID  U>  •  0(1) 

a 

<*•0 


t  -a  T  a  A  , 

4  •  v  0  x  •  x  y  ,  3  • 


ce",WW*  $'° . q*‘  • 


Then  it  follow*  laned  lately  fr^ss  (6.6)  that 

q-i  „ 

(6.7)  a  »  u  J  c  ,  a  •  0,...,q-l  . 

a  „  6 

9»0 


Hr  n.'w  prove  the  equivalence  of  our  five  condition*.  We  ahow  (i)  ■*  (ill*  (lii)  * 
(iv)  -  (l)  and  (l)  -  (v)  -  (Iv). 


(i)  -  (il)  Trivial. 


lii)  •*  (ill)  Since  1  la  A**  exist*.  1 1»  v*D  V*x  also  exist*,  a  •  0 . q-1  .  But 

a  *  •  a  »  *» 

v  u  >  0  •  and  hence  «  y'  •  0,  o  •  l,...,q-l  by  henna  (J.J). 

(Ill)  ••  (iv)  Consider  the  identity  (6.7).  since  the  Vanderaonde  matrix  q  l/‘(«  f>-), 
a.f  •  0,...,q-l  la  unitary  the  assumption  •  x  y  ’  *  0  .  o  *  l....,q-l  implies  that 
cQ  •  Cj  •  ...  •  cq.j  •  which  provea  (iv)  . 

(iv)  -  U)  If  (iv)  hold  then  c  -  C,  *  ...  »  c  and  (6.7)  iapliea  a.  -  ...  •  a  ,« 

0  1  q-i  1  q-1 

0.  Thia  eatabliahea  (1)  in  view  of  (6.5)  and  (6.6). 

(1)  —  (v)  Trivial,  since  vfx  >  0  and  u  >  0. 


(v)  ■*  (i)  let  »  -  ql  •  r  ,  0  ^  r  ^  q-1.  Then  (6.5)  iapliea 


ti*  Aq**r  «  -  x‘r’ 

•lT)  Z0.  x‘r*  Also 

.r  (0)  _ (r) 


r  •  0 . q-1 


r  (0)  (r)  „  ,  .q  (0)  (0) 

Ax  »  x  ,  r  •  0 . q-1  ,  A’x  -  x 

As  A1  is  a  direct  eua  nf  q  irreducible  and  prlaltive  aatricea  the  aaauaption  *  ^  0, 

x  *  Iapliea  that  11a  (A^)  x  •  x(  ’  ^  0.  Obviously  x 1  ’  ^  0  . 

t  — 

Mow  (v)  iapliea  that 

0  (1)  .  0  (0) 

*  Z  *  »  Ax  Z  x  , 

whence  x  ‘  »  x  and  thus  xr'  *  x  for  r  *  l,....q-l  .  So  1  in  A*x  ex10* 


and  (i)  follows 


In  what  follows  wv  give  necessary  and  sufficient  conditions  on  a  reducible  matrix 
A  to  satisfy  (e.2>.  To  do  so  we  need  a  few  more  graph  theoretical  concepts. 

Let  G  he  a  graph  on  <v>  •  (l,...,v).  Let  J  be  a  non-void  subset  of  (v>  . 

Then  i  ,  J  is  called  a  final  state  with  respect  to  J  If  for  any  f»  F  a  and  (a,f<)  r  G, 

ti  J  .  denoting  by  f(J)  the  set  of  all  final  states  with  respect  to  J  .  If  J  •  <v> 

then  i  is  called  a  final  state,  i.e.  (a,?)  <  G  implies  that  B  •  a.  Define 

d(g,J)  «  max'MS.a)  i  a  «  F(J1)  . 

If  J  •  <v>  then  write  d(g)  instead  of  d(B,<v>). 

let  A  _  be  a  reducible  matrix.  We  aasimr  that  A  is  In  the  rrobenlus  form  (4.1). 

As  in  14.  denote  by  G (A)  the  (reduced)  graph  of  A  .  Let  x  0,  x  f  0.  Partition 

x  conformably  with  A  given  by  (4.1).  That  is  x*  •  (x*^ . x*  ^).  The  support  of 

x  Is  the  set  supp  x  -  (aj , . . . ,  .i^  )  _  (1 . v!  such  that  x(j)  /  0  if  and  only  if 

i  *  »upp  x.  We  shall  always  assume  that  ,  1  *  l,...,s  have  been  listed  In  strictly 

ascending  order. 


(6.B)  Theorem 

Let  A  <  ,  f(A)  •  1.  Assume  that  A  Is  In  the  rrobenlus  form  (4.1).  Moreover 

if  Ai(  is  imprimltlvo  then  At i  Is  the  rrobenlus  form  (6.1)  let  «  >  0. 

*  *  3.  Then  (6.2)  holds  If  and  only  if  any  final  state  i  with  respect  to  the  support 
of  s  satisfies 


(I)  a  Is  s  singular  vertex  (l.e.  o(A  )  •  1) , 

on 

(II)  either  A  Is  primitive  or  A  and  x.  . 

-  s-i  —  aa  (a) 

(6.4)  . 

frocf ■  First  we  note  that 
(6.91  ( A_x)  • 

r 

Suppose  that  •>  .  F(sur?  x)  .  Then 


(a”x>  -  T  A(*'  X 

“  8  r  supp  X  38  (#1 


satisfy  the  condition  (iv)  of  1 


(Ax)  *  A  X,  . 

a  aa  ( a) 

By  the  definition  of  R(x)  and  r(x)  we  have 

r(AP'x)A*x  A**lx  B(AWx)AWx  . 
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So 

r (A**) A*  x,  ,  <  A**  *  X,  .  •  BlA^xlA*  X,  .  . 

aa  to)  —  aa  to)  -  aa  la) 

Hence,  since  A  irreducible,  by  (6.1), 
aa 

r(A*x)  «  r(A*  X,  ,)  «  0(A  >  <  B (A*  X,  .  )  *  B(A*x)  . 

—  aa  (a)  —  aa  —  aa  la)  — 

Assist#  now  that  (6.3)  holds.  Then  for  any  final  state  a  with  respect  to  supp  x  we 
Bust  have 


li»  r(A  x,  . ) 
aa  ( a) 


1  is  v  ( A  X  .  , ) 

aa  ( a) 


(s(A  )  •  1  . 
aa 


So  i  is  a  singular  vertex.  If  A  is  lapriaitive  then  the  condition  (v)  of  V 

aa 


6.4 


holds.  Hence  A  and  x,  ,  satisfy  (iv)  of  hesasa  6.4.  This  proves  one  direction  of 
aa  (a)  ' 

our  theorea. 

Ass use  now  that  if  a  <  Flsupp  x)  then  p(A  )  -  1  and  if  A  is  not  pnnltive 

aa  aa 

then  A  and  x,  .  satisfy  the  condition  llv)  of  I.eseaa  6.4. 
aa  (a) 

bet  1  s  «_  V.  Lot  d  -  d(fl,J).  By  our  assuertion  d^-1.  If  d  ■  —  ,  then 
.  If  d  ^0  ,  then 
-d 


(Ax) ^  •  0  ,  a-  1,3,... 


-<t  a 

•  U  x). 


a‘*'  X  »  0(1) 
6a  a 


where  %  •  (aiklf.a)  "  d*.  Clearly  B  c  Flsupp  x)  .  Thus,  to  show, 

(6.10)  li 

B  • 

it  is  enough  to  prove 

(6.11) 

for  a  #  F(supp  x),  X(8,.i)  -  d.  To  prove  (6.11),  let  0  be  the  Matrix  obtained  fro 


s  U  «l  >  0 


-d  . (a) 

1  *6a  *a'  0  • 


Y< 

-d.  (B) 

I  A.  X  •  B 

6a  a 


1  la  a 


l  o 


and  by 


(6.41 


p“0 


6a 


in 

all  other 

c%nt»n 

i  * 

0‘B*P’  AP 

x  . 

P“0 

6a  aa 

a 

f  r<-» 

8  to  a 

it  < 

P)  . 

0.  »  0 

6a 

1  in  A*"  x  -  v  >0. 

aa  a  a 

P  -  “ 
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(2.7)  that 


It  easily  follows  Item  lemma  (2.7)  that 

lla  a  A  *  "  T  U  V  »  0  . 

tvt  a  d  pi  j 

»  *  '  » 

Thus,  for  each  ,  1  ^  8  ^  v  ,  either  iA  K)^  *  0,  a  •  1,2,...  or  (6.10)  ia  satisfied. 

Pro*  this  (6.2)  follows  lamed  lately .  a 

We  remark  that  the  restriction  p(A)  «  1  in  Theorem  (6.8)  may  be  replaced  by  p(A)  >0. 

(6.12)  Ccrol lary i  bet  A  .  *”n  .  p(A)  -  1.  Assume  that  A  is  the  Frobenius  form  (4.1). 

bet  J  be  an  non-empty  set  of  (  v  > .  Then  for  any  x  ^  0  whose  support  is  the  set  J  , 

(6.2)  holds  if  and  only  if  for  all  final  state  ■  with  respect  to  J  ,  p(A  )  •  1  and 

30 

A  la  primitive.  • 

(6.U)  ■  orol  lary :  Lot  A  #  *  ,  c(A)  •  1  Aaauao  that  A  i»  in  the  Frobanlua  f on* 


Than  for  any  *  0  ,  *  t  t  (6.2)  hold*  if  and  only  if  for  aach  o,f( A  )  •  1 

*“  30 

and  A  ia  prlaitlva,  o  •  1 . .  * 

ao 


’*V  •  .• 


Aa  an  Application  of  our  results  we  <nvr  a  simile  proof  of  a  theorem  concerning  non- 


neqativs  solution*  y  of  (1  -  A)y  »  x  *or  aiven  *  _  0  .  for  l  a, 8  ^  v  we  shall 
say  that  has  access  to  a  In  C(A)  if  thete  la  a  (vath  fr<»  to  »  in  C  (A)  , 

wi*  A(8,o)  »_  -1. 

(7.1)  Thmjrcm  Let  A  I*'  '  with  p(A)  •  1  and  sunoose  that  A  is  in  the  froheniu* 
normal  form  (4.1).  Let  x  •  t  '  .  "hen  the  following  are  equivalent  i 
(1)  There  is  a  y  .  such  that  (I  -  A)y  •  *. 

(ii)  So  singular  vertex  p  has  access  in  CIA)  to  any  a  <  supp  x  . 

(ill)  l la  (I  ♦  ...  «  ab)x  exists, 

a  a  • 

(iv)  Ua  a  x  •  0  . 


Pur t  her . 

if  (ill)  holds 

and 

y  * 

lia 

a  •  • 

(I  ♦ 

A  ♦  ...  ♦ 

A*”)* 

,  then 

(X 

-  A)y  »  x  and 

(7.3) 

y8 

•  0  # 

if 

S 

does 

not  have 

jrc*n» 

to  any 

a 

•  supp  X  . 

(7.  J) 

ye 

>  o  • 

ir 

8 

has 

access  to 

a  f  supp 

X 

• 

Proof.  Let  S  **  •  I  ♦  A  *  . . .  ♦  A* 


_  If  1  '  (  ^  ,  then 

IS  a.  »  V  _(■> 

(7.4)  (S  x)„  •  )  S„  x  , 

P  *•  So  a 

1  >  supp  X 

and,  by  Corollary  (5.11),  for  A  -  klf.il  >_  -1, 

(7.5.1) 

a 

while  for  A(8.j)  •  -- 


,,  - (k*l)  (a) 

1  la  a  S  •  (1  »  0  . 

8a  So 


(7.5.11)  s‘"'  -  U-  -  O.a  -  1,2,1 . 

So  So 

Me  shall  prove  (1)  -  (U>  -  (iii)  -  (i),  (ill)  -  (lv)  -  (11)  . 

(1)  ■*  (11).  Suppose  that  (I  -  A)y  •  x  ,  where  y  __  0.  Then 
S<n'  x  •  (1  -  A**l)y  y  . 

Let  g  Sw  a  singular  vertex.  If  has  access  to  i  then  A  •  A(S.o)  ^  0  and 
by  (7.4)  and  (7.5) 


for  large  a 

(ii)  -  (iii) 


,  <*>  ,  1  0t*l)„ 

>•  >  (S  x)  „  >  -  a  U.  x 

8  -  S  —  3  So  a 

.  Hence  x  •  0  and  i  /  supo  x . 

Q 

Suppose  (ii)  holds  and  let  1  _  _  t. 


then  K 


exists,  for  a  ■  supp  * 


Since  A  S 


(tv)  *  (til.  Suppose  thst  (tv)  holds  but  thst  (it)  is  false.  Then  there  exists  a  slnqular 


8  and  an  >  .  supp  *  such  that  *<((,;»)  ^  0  .  Let  q  •  q(g,,,)  be  the  local  period  and 


•  A1  1.  Then  1  isi  B  x  -  0  ,  But  by  Theorem  (5.10)  for  all  suf 


flctently  larqe 


where  c  »  P  ,  and  x  f  0  .  This  is  a  contradiction,  and  the  Implication  is  proved 


To  complete  the  proof  of  the  theorem  observe  that,  for  y  *  1 1»  S 


n <  supp  x 

in  view  of  (11)  and  (7.5).  Since  U.  »  0  If  8  has  access  to  a  and  0. 

8a  8a 

otherwise,  w*  lassetiately  obtain  (7,2)  and  (7.J).  • 


The  equivalence  of  conditions  (l)  and  (ii)  in  Theorem  (7.1)  is  due  to  D.  M.  Carlson 


(1)  .  Ke  rrsnark  that  Carlson  also  showed  that  if  a  solution  y  of  (I  -  Aly  »  x  exists 


then  the  solution  satisfying  (7.2)  »nd  (7.1)  is  unique 
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Non-negative  matrices,  powers,  reduced  graph,  singular  distance 
iteration . 


Let  A  be  a  non-negative  n*n  matrix.  In  this  paper  we  study  the  growth 
of  the  powers  A1”  ,  m  »  1,2,3,...  .  These  powers  occur  naturally  in  the  itera 

tion  process 

x*-1*  «  AX(W>  ,  x<°’  >0 

which  is  important  in  applications  and  numerical  techniques.  Roughly  speaking, 
we  analyse  the  asymptotic  behavior  of  each  entry  of  A®  .  We  apply  our  main  re¬ 
sult  to  determine  necessary  and  sufficient  conditions  for  the  convergence  to  the 
spectral  r a diuso  f _ a _ of  certain  ratios  naturally  associated  with  the  iteration 
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